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Abstract 

Using special quasirandom structures (SQS's), we perform first-principles calculations studying 
the metastable bcc NiCu and FeCu alloys which occur in Fe-Cu-Ni alloy steels as precipitated 
second phase. The mixing enthalpies, density of state, and equilibrium lattice parameters of these 
alloys are reported. The results show that quasi-chemical approach and vegard rule can well 
predict the energetic and structural properties of FeCu alloys but fail to yield that of NiCu. The 
reason rests with the difference of bond energy variation with composition between NiCu and FeCu 
alloys induced by competition between ferromagnetic and paramagnetic state. Furthermore, the 
calculated results show that the energetic and structural properties of these alloys can well explain 
the local composition of the corresponding precipitates in ferrite steels. 

PACS numbers: 64.60.My, 61.43.Bn, 64.75. Op, 71.20.Be 
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I. INTRODUCTION 



The transitional metals and their alloys related to the magnetism attract much scientific 
interests. Enormous experimental and theoretical investigations were dedicated to deeper 
understanding of the nature of their properties which became complicated due to magnetism. 
Examples of these systems can be provided by artificial crystalline structures fabricated by 
film growth techniques and precipitated second phase in matrix. Specially, many impor- 
tant properties of materials, such as their mechanical strength, toughness, creep, corrosion 
resistance, and magnetic properties are essentially controlled by precipitated particles of a 
second phase. Therefore, understanding of these alloys is also desirable from the views of 
applications. 

The Cu-rich precipitates, which are commonly found in alloy steels, have good strength- 
ening effect on steels. The high-strength low-alloy steels strengthened by copper rich pre- 
cipitates also retain the impact toughness, corrosion resistance, and welding properties[l-6|. 
However, it is also confirmed that the presence of Cu-rich precipitates in reactor pressure 



vessel(RPV) steels is the origin of embrittlement [7[, which limit the reactor operating life. 
The embrittlement effect can be enhanced by the content of Cu or Ni of alloyed steels con- 
taining of both Cu and Nijs-JJ]. Thus, the mechanism of strengthening and embrittlement 



of Cu precipitates attract much interest 



13|-|22||. The structure and composition of precip- 



itated phase also became important, and it had been studied by many experiments, such 
as atom probe field ion microscopy(APFIM) [23-25], small-angle neutron scatterin gfSA NS) 



26| . high- resolution and conventional electron microscopy(HREM and CTEM) 27|, |28|, etc. 



It was observed that small Cu-rich precipitates with diameters less than about 5 nm have 
a meta-stable body-centered cubic (bcc) structure and are coherent with the a-Fe matrix 



26|-|28| in the initial stage of segregation. In addition, it was also found that Ni also occurs 



in the Cu-rich precipitates beside Cu and Fe. For HSLA steels, Ni was observed at the co 



herent matrix/Cu-precipitate heterophase interfaces [29|, l30|, l36|, l37| . For RPV steels, it was 



also observed Ni appeared in the Cu-rich precipitates at very initial stage and was rejected 



from the core with the precipitates growth 
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The bcc FeCu alloy has been studied intensively by theoretical investigations for well 
understanding of precipitated phase. First-principles calculations based on the cluster ex- 
pansion framework gave the composition range of mechanical stability of FeCu random 
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alloy 351]. A thermodynamic equilibrium analysis was used to investigate the composition 
dependence of Gibbs energy of FeCu alloy as well as the interfacial energy between precipi- 
tates and matrix 361] . In the meanwhile, the vibrational energies of dilute FeCu alloys were 



investigated by first-principles calculations in combination with thermodynamical model- 
ing, which show that the vibrational energies can stabilize the alloys 37|. The dependence 
of magnetism on the structural characteristics of Cu nucleation and the electronic struc- 
ture for bcc Fe2;Cui_2;(x > 0.75) systems were also well investigated |38f]. However few 
systematic investigation of bcc NiCu alloys was reported. The bcc NiCu, like bcc FeCu, is 
not only a typical system related to magnetism in solids, but also became interesting be- 
cause of its important contribution to strengthening and embrittlement of the ferrite steels. 
Therefore, systematic investigations of their properties at electronic level have both founda- 
tional and engineering significance. In this paper, the composition-dependent bcc Ni^Cui-^; 
and Fea;Cui_a; random alloys are investigated by employing special quasi-random structures 
(SQS's) in the frame of first-principles calculations. 



II. MODELING AND THEORETICAL METHODS 



The SQS method was 



of mean-field theories 



41 



proposed by Zunger and Wei et al to overcome the limitations 



421 ] ■ without the prohibitive computational cost associated with 



directly constructing large supper cells with random occupy of atoms. The SQS method 
was extensively used to study the properties of semiconductor alloys which are all fcc-based 
system. Recently, this method was also used into both fee and bcc transition metal systems 
and was proved to be useful 



43 



44| . For a binary substitutional alloy, many properties are 



dependent on its configuration. A binary AB substitutional alloy with a lattice of sites 
has 2^ possible atomic arrangements, denoted as configurations a. The measurable property 
< E > represents an ensemble average over all 2^ configuration a, < E >= p{cr)E{a). 
The structure a can be discretized into its component figure /, which is characterized by a set 
of correlation functions 11^ ,^. Therefore, < E > can be rewrite as < E >= ^ 11^ ,„£(/)• The 
SQS's are special designed A^ atoms periodic structures whose distinct correlation functions 
IlA;,m best match ensemble averaged (11^^^) of random alloys. 



The choice of SQS's are critical for the calculations 



45| . In this paper, SQS's containing 



32 atoms are constructed for bcc alloys, which are shown in Fig. 1. The vectors of lattice are 
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ai = (1.0,-2.0,0.0)ao,4 = (0.0, -4.0, 2.0)ao, a| = (-2.0, 0.0, -2.0)ao, respectively, and Oq 
is the lattice parameter of the bcc unit cell. The occupations of these sites for the SQS's at 
X = 0.25 and x = 0.5 are given in Table 1, and their structural correlations function 11^^^ 
compared with ideal random alloy correlation functions are given in Table 2. As can be 
seen, the quality of the SQS's used in these calculation is reasonably good. 

First-principles calculations are performed using the density functional theory 
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48| as 



implemented in the Vienna ab initio simulation package(VASP) 49|]. The generalized gra- 
dient approximation(GGA) 50|, l5l(| is used for the exchange correlation functional. The 
interaction between core and valence electrons is described with the projector augmented 
wave (PAW) potential 52i l53| . The equilibrium structures were determined, up to a preci- 
sion of 10~^ eV in total-energy difference and with a criterion that required the force on each 
atom to be less than 0.01 eV/A in atomic forces. The convergence tests with k-point and 
energy cutoff are presented in table 3. The k-point meshes 6 x 2 x 4 for the Brillouin zone 
integration are used for SQS's, which is equivalent k-point sampling with Monkhorst-Pack 
mesh 11x11x11 for bcc primitive cell. An energy cutoff of 280 eV is applied in all cases. 



III. RESULTS AND DISCUSSIONS 



Firstly, basic properties of bcc Fe, Ni, Cu are presented. The lattice constant and the 
magnetic moment for bcc Fe at the equilibrium volume are 2.83 A and 2.16 //fi/atom, 
respectively. The lattice constant and bulk modulus for bcc Cu at the equilibrium volume 



are 2.89 A and 130 GPa. These results are in agreement with earlier studies 35l. |39|. The 
magnetic moment for bcc Ni at the equilibrium structure with lattice of 2.81 A is 0.55 
/ifi/atom. The calculated result is we 



predicted bcc Ni is ferromagnetic 54 



1 consistent with previous works using GGA, which 



55| . Previous theoretical works using local density 



approximation (LDA) predicted bcc Ni is paramagnetic 39|. However, it is confirmed that 



bcc Ni do possesses a magnetic moment of 0.52 /i^/atom by experiments more recently |56[|. 
In addition, it was proved that the GGA calculations can describe satisfactorily lattice and 
elastic constants of Ni in their observed structures jssl. 

The mixing enthalpies(Ai^f) of bcc NiCu and FeCu random alloys obtained by first- 
principles SQS method are shown in Fig .2. The mixing enthalpy is heat quantity ab- 
sorbed(or evolved) during mixing two elements to make homogeneous solid solution, which 
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is defined as: 

AH = Eab - xEa - (1 - x)Eb (1) 

where E^, Eb, and E^b are the total energy of A, B, and AB alloy, and x is the compo- 
sition of alloys. In the frame of quasi-chemical approach[57], the mixing enthalpy is only 
determined by the composition- independent bond energy between adjacent atoms. Hence, 
mixing enthalpy can be written as: 

AH = nx{l-x), (2) 

where fl = Naz{eA-B — \{,£a-a + ^b-b)) is interaction parameter, is Avogadro's number, 
e is bond energy, and z is the number of bonds per atom. Thus, the composition-independent 
bond energy between atoms would result a parabolic variation for mixing enthalpies with 
compositions. It can been seen from Fig. 2, the variation of mixing enthalpies of FeCu 
with composition is parabolic, while that of mixing enthalpies of NiCu exhibit a strong 
asymmetry and change sign as function of concentration. These effect were also found for 
random fee NiCu alloys in previous calculations |4o|. 

To study the effect of degree of disorder on the mixing enthalpies of NiCu, we calculate 
the ordered structures comparing disordered structures. Three 4-atoms-super-cells are used 
to simulate order alloys with composition x =0.25, 0.50, 0.75. The mixing enthalpies of these 
order alloys are represented by solid squares in Fig. 2(a). It shows the mixing enthalpies of 
this ordered structures of NiCu are all lower than those of disordered alloy with the same 
extent, and the trend of mixing enthalpies variations with compositions are the same as that 
of disorder alloys. The mixing enthalpies of order alloy also exhibit a strong asymmetry. 
Therefore, the degree of disorder is not the reason of the asymmetry of mixing enthalpies 
for NiCu alloys. 

We turn to investigate bonding property of NiCu and FeCu from electronic structures 
for understanding strong asymmetry variation of mixing enthalpies with composition. The 
bond energy is critical factor that determines the mixing enthalpies. The variation of the 
bond energy can be reveal from comparison of the density of states(DOS) of Fe^Cui-i; and 
Nij^Cui-a;. Since the bond energies of transitional metal are determined by the coupling of 
d-band, the d-bands of Cu, Fe, Ni in FcajCui-x and Ni^^Cui-j^ with different compositions are 
plotted in Fig 3. One can easily understand the difference of mixing enthalpies magnitude 
between FeCu and NiCu from electronic structure. In comparison with the occupation of 
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Ni, the occupation of Fe in the alloy is changed significantly from that of pure Fe(see DOS 
around 2 eV in Fig. 3). Therefore, the energy of Fe-Cu bond is much larger than that 
of Ni-Cu bond, and the mixing enthalpies of CuFe are much larger than those of CuNi. 
Furthermore, bond energies in alloys vary with compositions can be also refiected by the 
DOS. The DOS peaks of Ni vary with composition, while that of Fe in alloys with different 
compositions are the same. These findings indicate that bond energy e of Ni-Cu varies 
with composition. It can be known from equation (2), the variation of bond energy induce 
composition-dependent interaction parameter fi, which result that the variation of mixing 
enthalpy with composition is not parabolic. These are consistent with calculated results in 
Fig 2(a). Hence, the complex variation trend of mixing enthalpies of NiCu derive from the 
composition-dependent bond energy. 

In addition, one can infer that the magnetic moment of Ni atom varies with the com- 
position and the magnetic moment of Nio.25Cuo.75 disappears. A more clear presentation of 
magnetic moment of alloys with different composition is given in Fig. 4. These phenomena 
can be well understand by using stoner criterion which succeeds in explaining the magnetic 
property of transitional metal alloys dictated by the filling of the d-band39l|. The Stoner 
criterion states that ferromagnetism appears when the gain in exchange energy is larger 
than the loss in kinetic energy. Therefore, there is always a competition between ferro- 
magnetic(FM) and paramagnetic(PM) solutions, and magnetic properties are determined 
by the state which has lowest energy. The DOS at the Fermi level is a key factor to affect 
the competition between FM and PM state. A larger peak of DOS at Fermi level induces a 
larger exchange energy and prefers to split into FM state. Fig. 5 illustrates the unpolarized 
DOS of NiCu and FeCu alloys. It is shown in Fig. 4 and Fig. 5, the trend of magnetic 
moment of the system well consistent with the variation of DOS at Fermi level. 

For these alloys, the d electrons increase with the Cu concentration, which result the 
Fermi level of Ni(Fe) be promoted. The key point is, as shown in Fig. 5, since the DOS 
at Fermi level of NiCu is much smaller than that of FeCu, the Fermi level of NiCu change 
much rapidly with Cu concentration than that of FeCu which induce that the peak of DOS 
of NiCu is much sensitive to the alloy composition than that FeCu. Hence, the magnetic 
moment of Ni change with the variation of DOS at Fermi level when the Cu concentration 
is changed, and it disappears when the DOS at the Fermi level is 0. Now, we can conclude 
that the strong asymmetry of mixing enthalpies of NiCu with composition is induced by 
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band shift because of the competition between FM and PM state. 

It is found that the composition-dependent electronic interactions of NiCu alloys also 
affect the structural property. As shown in Fig. 6, the calculated volume per atom of FeCu 
alloy well obeys the vegard law, while that of NiCu deviates from what the vagard law 
predicts. According to the vegard law|58|, the lattice constant has a linear relationship with 
its composition, i.e. a{AxBi^x) = xa{A) + (1 — x)a{B). This implies that the atom can be 
considered as rigid body whose volume is unchanged with its chemical environment. How- 
ever, the composition-dependent electronic interaction in NiCu result the atomic volumes 
depend its composition, which induces bowing effect. 

As shown in Fig. 7, the partial atomic pair correlation functions give more detailed 
information of the alloy structures. It can show that how large is the degree of distortion of 
equilibrium alloy structure from the perfect bcc structure. For bcc structure, the numbers 
of 1st, 2nd, 3rd, 4th, 5th neighbors are 8, 6, 12, 24, 8. As is shown in Fig. 7, both NiCu and 
FeCu alloys are still bcc structure and the lattice just has a minor change, but the NiCu 
alloy is more close to a perfect bcc structure. Since the nominal number of minority d-band 
holes of Fe is more than that of Ni by 2 electrons, the partial charge effect of Fe from Cu 
is stronger than that of Ni. The bond length can be affected by partial charge. Hence, the 
distance between atomic neighbors of FeCu more sensitive to their chemical environment 
than NiCu, which induce a larger distortion of FeCu from bcc. 

The strain effect is also very important for the precipitated phase in ferrite matrix. We 
have calculated the strain energy with lattice constant from 2.80 A to 2.90 A as shown in 
Fig. 8, which covers the lattice constant range of bcc Fe and Cu. It is shown that NiCu 
alloys have smaller strain energy at lattice of bcc Fe, while FeCu alloys have smaller strain 
energy at lattice of bcc Cu. At a lattice constant of Fe matrix with 2.83 A, the strain energies 
are 6.52, 0.78, 0.51 meV/atom for Ni^Cui_^(x=0.25, 0.5, 0.75) alloys, and 20.71, 13.51, 7.13 
meV/atom for Fei.Cui_a;(x=0.25, 0.5, 0.75). The strain energy of bcc Cu in ferrite matrix is 
21.6 meV/atom. These indicate that Ni and Fe can lower the bcc Cu strain energy in ferrite 
matrix, while the Ni is better that Fe in role of reducing strain energy. 

From the energetic and structural properties of alloys, the precipitated phase can be well 
understood. For HSLA steels, it was observed that NiCu alloys at the matrix/precipitates 
heterophase interfaces area|29i]. The calculated results show that the mixing enthalpies of 
FeCu are positive, which indicate that the FeCu alloys are unstable and Cu atoms prefer 
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to segregate to form precipitates. The mixing enthalpies of NiCu alloys are negative when 
Ni concentrations below 30 at. %, indicating that NiCu alloys with lower Ni concentration 
are energetic favorable. The element of Ni appears in the matrix/precipitates heterophase 
interfaces area, which is mainly induced the factors relate to strain: i) The volume per 
atom of NiCu is similar to that of Fe matrix rather than pure Cu, even rather than that 
of FeCu alloy, and the stain energy can be much lower by Ni. ii) The structure of NiCu 
is more close to the bcc structure, which induce little mismatch of lattice. These indicate 
that NiCu alloys in the interface of matrix/precipitate can contribute lower interface energy 
from reducing shear stress. However, the experimental peak concentration of Ni is about 3 
at. % in the Cu-rich precipitates which is lower than the Ni concentration of most favorable 
NiCu alloy with concentration of 20 at. %. This indicates that the distribution of elements 
distribution may be affected by the size effects and vibrational entropy in addition. In RPV 
steels, the Ni in Cu-precipitates be rejected from the core after thermal-aging or neutron- 



irradiated 
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34j . and the experimental peak concentration of Ni is among 15-20 at. %, 



which is well consistent the composition of energetic favorable alloys more closed. These 
indicate that the element distribution of bcc precipitated phase at very initial of segregation 
stage can be well understood from the energetic and structural properties. 



IV. CONCLUSION 



In summary, we propose three 32-atom SQS supercells to mimic the pair and multisite 
correlation functions of random CuFe and NiCu bcc substitutional alloys which occur in 
Fe-Cu-Ni alloy steels as precipitated second phase. Those SQS's are used to calculate the 
mixing enthalpies, density of state, and lattice parameters of the metastable random alloys. 
The results show that quasi-chemical approach and vegard rule can satisfactorily predict the 
mixing enthalpies and structure parameters of FeCu alloys but fail to accurately yield that of 
NiCu. As can be obtained from the analysis of electronic structure, the magnetism induced 
bond energy variation with composition is the reason that quasi-chemical approach and 
vegard rule fail to predict the properties of NiCu alloys. Furthermore, the calculated results 
can well explain the previous experimental observation of local composition of coherent 
Copper-rich precipitates containing Nickel and confirm that segregation of Ni is drove by 
thermaldynamic and chemical factor. These suggest that the properties of random alloys 



8 



have important implications to better understand the multi-component precipitates. 

In this work, we present the intrinsic bulk properties of metastable bcc FeCu and NiCu 
random alloys to understand the structure of Cu-rich precipitates. Since Cu-rich precipitates 
have important roles on properties of alloy steels, understanding the formation mechanism 
is highly desirable for altering the properties of steels by controlling Cu-rich precipitates. 
Further investigation would simulate realistic interface to calculate the interfacial energy 
of precipitates/matrix and the diffusion properties of Cu and Fe atom affected by Ni shell, 
which will clarify the role of Nickel on the evolution of the Cu-rich precipitates. 
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TABLE I: Atomic coordinates and occupation of the 32 atoms bcc SQS. 



Atomic coordinates 



AB 

0.0-2.0 0.0 
0.0-1.0 0.0 
-1.5 -2.5 -0.5 
-0.5-3.5 0.5 
-0.5-5.5 0.5 
-1.5 -0.5 -1.5 
-0.5 -2.5 -0.5 
0.0-4.0 1.0 



-0.5 -0.5 -0.5 
-1.0 -3.0 -1.0 
-1.5 -3.5 -0.5 
0.0-3.0 1.0 
0.0-5.0 1.0 
-1.5 -1.5 -1.5 
-0.5 -3.5 -0.5 
-1.0-3.0 0.0 



A: 0.0-4.0 0.0: 
A: -0.5 -1.5 -0.5: 
A: 0.5-5.5 1.5: 
A: -1.0 -6.0 0.0: 
B: -0.5 -4.5 0.5: 
B: 0.5-3.5 0.5: 
B: -1.0 -2.0 -1.0: 
B: 0.0-6.0 1.0: 



A: 0.0-2.0 0.0: 
A: -0.5 -3.5 -0.5: 
B: 0.0-4.0 0.0: 
B: -0.5 -0.5 -0.5: 
B: -1.0 -3.0 -1.0: 
B: 0.5-5.5 1.5: 
B: -1.0 -6.0 0.0: 
B: -1.0 -5.0 0.0: 



AB, 

0.5-2.5 0.5; 

0.0-3.0 0.0; 
-1.5 -0.5 -1.5; 
-1.0 -4.0 -1.0; 
-0.5 -2.5 -0.5; 
-1.5 -3.5 -0.5; 
-1.0-4.0 0.0; 
-0.5-4.5 0.5; 



0.5-3.5 0.5: 
0.0-4.0 1.0: 
-1.5 -1.5 -1.5: 
-1.0 -2.0 -1.0: 
0.0-6.0 1.0: 
-1.0 -3.0 0.0: 
-0.5-3.5 0.5: 
-0.5-5.5 0.5: 



-1.0 -4.0 -1.0 
0.5-4.5 1.5 

-0.5 -2.5 0.5 

-1.0-5.0 0.0 
0.5-2.5 0.5 

-1.0 -1.0 -1.0 
0.0-3.0 0.0 

-1.0-4.0 0.0 



0.0-1.0 0.0 
-1.5 -2.5 -0.5 
-1.0 -1.0 -1.0 
-0.5 -1.5 -0.5 

0.5-4.5 1.5 
-0.5 -2.5 0.5 

0.0-3.0 1.0 

0.0-5.0 1.0 
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TABLE II: Comparison of correlation function Ilk,m of 32 atoms bcc SQS and ideal value of random 
alloy. 





AB 




AB3 






random 


SQS 


random 


SQS 


Ho 1 








0.25 


0.25 




U 


U 


U.ZO 


U.ZO 


n2,3 








0.25 


0.25 


n2,4 








0.25 


0.25 


n2,5 





0.17 


0.25 


0.04 


n2,6 








0.25 


0.25 


n3,2 





0.03 


0.13 


0.13 


n4,2 








0.06 


0.08 



TABLE III: Total energies of bcc primitive cell of Fe, Ni and Cu as a function of the number of k 
points and the cutoff energy. 



Monkhorst-Pack mesh 


energy cutoff (eV) 


bcc Fe 


bcc Cu 


bcc Ni 


3x3x3 


280 


-8.15 


-4.13 


-5.27 


5x5x5 


280 


-8.12 


-3.71 


-5.33 


7x7x7 


280 


-8.17 


-3.70 


-5.39 


9x9x9 


280 


-8.17 


-3.70 


-5.37 


11 x 11 x 11 


280 


-8.17 


-3.70 


-5.37 


13 X 13 X 13 


280 


-8.17 


-3.70 


-5.37 


11 X 11 X 11 


180 


-6.99 


-1.97 


-4.05 


11 X 11 X 11 


260 


-8.15 


-3.68 


-5.36 


11 X 11 X 11 


300 


-8.17 


-3.70 


-5.37 
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(a) 



(b) 



FIG. 1: Crystal structures of the SQS-32 supercell in their unrelaxed forms. Blue and red spheres 
represent A and B atoms, respectively, (a) AB (b) AB3 
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FIG. 2: Mixing enthalpy of alloys: (a) disordered and ordered Nia;Cui_a;, (b) disordered Fei.Cui_ 
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FIG. 3: Comparison of the d-band between pure metals and alloys. The d-band of Ni, Fe, Cu in 
pure Ni, Fe, Cu, NiCu alloy, FeCu alloy are plot in each panel. 
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FIG. 4: Magnetic moments per atom in NiCu and FeCu alloys. 




FIG. 5: The unpolarized calculated density of state(DOS) of NiCu and FeCu alloys. The solid, 
dash, dot line denote for Ni(Fe)a;Cui_a;, x = 0.25,0.50,0.75, respectively. 



16 



o 

< 

s 

> 



12.2 
12.0 
11.8 
11.6 
11.4 
11.2 
11.0 



1 


■ Fe Cu 

X 1-X 




• Ni Cu 

X 1-X 















0.00 0.25 0.50 0.75 
Atomic Fraction 



1.00 



FIG. 6: Atomic volume versus composition in bcc Ni^.Cui_2: and Fe^^Cui-^^ alloys. The circles and 
squares denote the calculated results of SQS method. The solid lines denote the results predicted 
by vegard rule. 
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FIG. 7: Partial atomic pair correlation function gij{r) in NiCu and FeCu alloys. 
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FIG. 8: The strain energy of NiCu and FeCu as a function of lattice constant. 
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